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Abstract 

We consider an integrable model which describes light beams prop- 
agating in nonlocal nonlinear media of Cole-Cole type. The model 
is derived as high frequency limit of both Maxwell equations and the 
nonlocal nonlinear Schrodinger equation. We demonstrate that for a 
general form of nonlinearity there exist selfguided light beams. In high 
frequency limit nonlocal perturbations can be seen as a class of phase 
deformation along one direction. We study in detail nonlocal perturba- 
tions described by the dispersionless Veselov-Novikov (dVN) hierarchy. 
The dVN hierarchy is analyzed by the reduction method based on sym- 
metry constraints and by the quasiclassical 9— dressing method. Qua- 
siclassical 9— dressing method reveals a connection between nonlocal 
nonlinear geometric optics and the theory of quasiconformal mappings 
of the plane. 

Key words: Nonlocal Nonlinear Optics, Cole-Cole Media, Phase 
Singularities, Dispersionless Systems, Quasiconformal Mappings. 



1 Introduction. 

The optics studies phenomena of the propagation of the electromagnetic 
waves through a dielectric medium in absence of currents and charges pp. 
In such the Maxwell equations assume the form 

dU 

VAH- — =0 (1.1a) 

dt ^ ' 

VAE + — = (1.1b) 
at 

V-D = (1.1c) 

V-B = 0. (l.ld) 

where x, y, z are the spatial coordinates, V = (dx-, dy, dz) is the gradient and 
t is the time. For sake of simplicity we have set the light speed c = 1. The 
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vectors E and B are the electric and magnetic fields respectively, while the 
displacement vector D and the magnetic induction H contain the informa- 
tion about the response of the medium when an external electromagnetic 
field is applied. In general, D and H are certain functions of E and B and 
they are specified by the so-called constitutive relations. In several physically 
meaningful cases they can be written down as follows 

D = eE B = iiU, (1.2) 

where electric permittivity e and the magnetic permeability fi depend on the 
coordinates and on the fields. In the case in which s and fi depend at most 
on the coordinates, the Maxwell equations are linear and provides us with an 
exhaustive description of a very broad class of phenomena in physics ^EUHj. 

The nonlinear optics deals with the class of media such that e and /i 
depend on the fields. In this case, the Maxwell equations are nonlinear and 
construction of their solutions is, in general, a very challenging problem. In 
the following, we will focus only on the electric part of the field, assuming 
the magnetic response to be negligible, i.e. fj, = 1. This condition is realized 
in many experimental situations. 

By definition, a dielectric medium is referred to as nonlinear if the di- 
electric function e is a certain function of the electric field e = e (E). 

The solution of the Maxwell equations for an arbitrary form of the func- 
tion e is, of course, a hard problem. Nevertheless, several very interesting 
cases, such as quadratic and cubic nonlinearities, are amenable by exact 
methods. Moreover, they are connected with a relevant phenomenology 
such as the higher harmonic generation (quadratic nonlinearity) and the 
soliton production (cubic nonlinearity) 1^. 

In the present paper we consider a model for light beams propagation 
in the limit of high frequency. This model can be derived both as high fre- 
quency limit of the Maxwell equations and from the nonlinear Schrodinger 
(NLS) equation which describes paraxial light beams. We will show that it 
is reasonable to assume the nonlocality to be weak because of the quickly 
oscillating fields. In particular, we can separate the pure nonlinear con- 
tribution from the higher orders nonlocal perturbations. Leading order is 
represented by the standard eikonal equation. Note that in the standard 
theory higher order terms contain the amplitude of the electric field as well 
as the phase S and consequently in virtue of nonlocality and nonlinearity 
they are rather complicated. Nevertheless, there exist a type of media where 
it is possible to consider a set of nonlocal perturbations along one direction, 
say z, which are not mixed up with the wave contributions. These are the 
Cole-Cole media. They are characterized by the Cole-Cole dispersion law 
for the dielectric function of the form e = eq + 0(a;~") for < a < 1 and 
cj ^ oo. We would like to stress that a large variety of solid and liquid 
polar media obeys the Cole-Cole dispersion law j^. In this case we perform 
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an asymptotic expansion of all fields with respect to the small parameter 
a;~° separating first nonlocal correction to the phase S from the higher or- 
der wave corrections. In our derivations a suitable slow dependence along 
z— direction is assumed. Nonlocal corrections are described by polynomials 
in Sx and Sy. In particular, a one-to-one correspondence among nonlocality 
and polynomials degrees is realized. 

In the a; — > oo limit separation between the nonlinear response and the 
nonlocal effects along z suggest to study first the properties of the light beam 
on the xy— plane (transverse equations) and then the "evolution" along z. 
Transverse equations are the two-dimensional eikonal equation and the ge- 
ometric optics limit of the Poynting vector conservation law. It is shown 
that for a general class of nonlinear responses obeying the so-called ellip- 
ticity condition there exists self-guided light beams. Moreover, several light 
beam features are determined in terms of the properties of quasiconformal 
mappings on the plane via the Beltrami equation. Propagation of helicoidal 
wavefronts and its first degree nonlocal perturbations are also discussed. 
Thereafter, we focus on the special class of perturbations which preserve the 
phase inversion symmetry S — > —S similar to the eikonal equation. Such 
perturbations are shown to be described by an infinite set of integrable non- 
linear partial differential equations (PDEs). It is the dispersionless Veselov- 
Novikov (dVN) hierarchy. First non-trivial equation of the hierarchy can 
be formally obtained by the slow variable expansion of the Veselov-Novikov 
equation, which has been introduced in as {2 + l)D integrable generaliza- 
tion of Korteweg-de-Vries equation. The dVN hierarchy is of interest since 
it is treatable by different approaches. More specifically, we discuss the re- 
duction method based on the symmetry constraints. It allows us to reduce 
the dVN equation to a 1 -|- 1— dimensional hydrodynamic-type system and 
it is effective for construction of explicit solutions. Moreover, we study the 
eikonal equation and the dVN hierarchy using the quasiclassical 9— dressing 
method. It provides us with a general approach to construct and analyze 
the whole hierarchy. 

It is worth to note that the quasiclassical 5— dressing method establishes 
a remarkable connection between the nonlocal nonlinear geometric optics 
and the theory of quasiconformal mappings of the plane. 

The paper is organized as follows. The model under consideration is de- 
rived from the Maxwell equations in Section [21 and starting with NLS equa- 
tion in the Section |3J In the Section |1] we analyze the transverse equations 
in connection with the Beltrami equation and minimal surfaces equation 
which gives the helicoidal wavefronts. Nonlocal perturbations described 
by the dVN hierarchy are derived in the Section [5J After a short review 
about integrable system in the Section IHl we discuss the hydrodynamic-type 
reductions obtained using the symmetry constraints in Section [7| and the 
application of the quasiclassical 5— dressing method to the eikonal equation 
and to the dVN hierarchy in Section |H1 Some concluding remarks close the 
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paper. 



2 The general model. 

Let us consider a medium with the magnetic permeabihty /i = 1. In such a 
case the Maxweh equations p.lj) imply that 

VAVAE+— ^=0 (2.1) 

For time oscihating solutions, of the form 

B{x,y,z,t) = -E{x,y,z) e'^' 
Y>ix,y,z,t) = -D{x,y,z) e'"*, 

equation (|2.1j) looks like 

VAVAE-cj2d = 0, (2.2) 

or, equivalently, 

V^E + Lo^-D - V {V ■ E) = 0. (2.3) 

Once the constitutive relation for D is assigned equation (|2.3|1 determines 
the electric component of the field. 

The constitutive relation depends on the physical properties of the medium. 
Let us assume that the displacement vector D can be splitted in a local part 
D/, and a nonlocal one D/vr as follows 



where 



D = Di + D^ (2.4) 



= e (/(x)) E(x) (2.5) 



Djv = / (x' - X, a{uj)) N (/(x')) E(x') dx' . 



(2.6) 



The quantity / = |Eo| is the intensity of the electric field, e{I) and (/) 
are, respectively, the local and nonlocal nonlinear responses and the function 
i2(x— x') is the nonlocal distribution. The parameter a{uj) governs the width 
of the nonlocal response in different regimes of the frequency. In particular, 
in the high frequency limit it is reasonable to assume that the external field 
is not resonant with the proper oscillations of the particles of the medium 
in such a way that the response becomes local, i.e. 

lim R (x' - X, a{Lo)) = (5 (x' - x) , (2.7) 
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where 6{x.' — x) denotes the Dirac (5— function. In order to explain this 
behavior, one can consider a naive model where the single particle of the 
medium is described by a one dimensional forced oscillator 

X + ujqx = /o sin {cat) . (2-8) 

As usual, the dot denotes the total time derivative, cuq is the proper oscilla- 
tion frequency of the particle. The left hand side contains the forcing term 
represented by a linearly polarized electric field oscillating with frequency 
u). The solution of equation (|2.8I) is of the form 

fo 

c\ sin (wo*) + C2 cos (wq^) H 5 ^ sin [uot) uj^ ^ uj 

ci sin {ojQt) + C2 cos {ojQt) t cos {oj^t) ujq = uj 

We remind that we are interested in the non-resonant regime loq ^ u). In this 
case, it is easy to see that for a; —> 00 the forcing contribution disappears 
and the particle behaves as the harmonic oscillator. In this regime one 
expects that the effect of the external field, up to higher corrections, does 
not propagate far from the point considered and the response tends to be 
as local as much a; is larger. 

Substituting the expression 1)2. 4|) into equation ()l.lc() one gets 

V-E = -V(log£)-E- ^ _ (2.10) 
Thus, equation (|2.3j) takes the form 

V2E + a;2(eE + DAr) - V ^-Vloge-E - ^ '^^ ^ =0. (2.11) 

Let us introduce a general model of weak nonlocality. It is given by a 
nonlocal distribution function of the following form 

(x - x') = pi^n (x) (x - x') (2.12) 

where the sum on the repeated indices is assumed. The distributions S^^'"^'""^ 
are the Dirac 5— function derivatives. They are standardly defined via 

Jm.3 ^ ) J \ ) y > dx[dxl^dxf ^ ' 

where we set x = {x\^X2^x-i). 

Using the weak nonlocal distribution (|2.12|) . one can rearrange the ex- 
pression on the nonlocal contribution Djv as follows 

D;v = r(0)(x)7V (/(x)) E(x) + rf^ (x) d.,, {N (/(x)) E(x)) 

+ rl'2d.,d,^ {N (/(x)) E(x)) + rl'2n (x) d,,d,^d,„ {N (/(x)) E(x)) + . . . 

(2.14) 
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In principle the coefficients r^") might depend on the frequency. In particu- 
lar, we assume the following power dependence 



r^'^) = . (2.15) 



We consider a model where the function N depends on the intensity accord- 
ing to the formula 

iV(/(x)) = ii^ (2.16) 

where a is a real positive constant parameter. 

Under the assumptions mentioned above we perform the geometric op- 
tics (semi- classical) limit of equation p.ll|) . As usual, let us represent the 
electric field in term of the phase S 

E(x) = Eo(x) e^-^^W. (2.17) 

In high frequency limit, uj~" is the small parameter with respect to which 
we may consider the following asymptotic expansions 

I = Io + uj-^h + O (a;-2") (2.18) 
e (7(x)) = 80 (/o(x)) + LJ-^ei (/o(x), /i(x)) + O {u;-^^) . (2.19) 

Evaluating Djv in the limit of high frequency, one gets 

Djvf«eo(/o(x))Eo(x)Pe^'^^W, u; ^ oo (2.20) 

where 

V = f^°^ -I- if I ■* Sxi — SxiSx^ — ifl2n^XiSxmSxn + • • • , 

l,m,n = 1, 2, 3. 

Using the expression 1)2. 2U() . we get 

V -Bn ^iioV So (/o(x)) Eo(x) • e^'^^W (2.21a) 
V (V • D^) « -a;2 V eq (/o(x)) (Eq • V^) VS e^'^^W. (2.21b) 



Recall that the vector VS is perpendicular to the wavefront and provides us 
with the light rays direction. Imposing the so-called transversality condition 

Eo-VS = (2.22) 

we restrict ourselves to the solutions such that the electric field E is perpen- 
dicular to the light rays. In this case the terms (|2.21a|) and (|2.21b|) do not 
contribute to the high frequency limit of the equation ()2.11|1 . 
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Moreover, we look for z— variable slowly depending solutions of the following 
form 

S = kz + S{x,y,T) (2.23) 

where r := z/ijO°^ is the "slow variable". The high frequency limit of the 
equation (|2.11|) at the leading order gives 

Sl + Sl = 4u (/o(x)) (2.24) 

where 4u (/o(x)) = ffo (-^o(x)) — k"^ ■ 

If the parameter a is such that < a < 1 we get an intermediate con- 
tribution between the pure geometric optics order and the first correction 
containing the amplitude of the electric field. In order to calculate it, we 
note that 

Sl = {k + oj-'^Srf = k^ + 2kuj-''Sr + oj-^'^Sl (2.25) 

where we kept into account that 

d d 

The a;^~° order term in equation (|2.1H) is 

2kSr + V (x, y, r, 5^, Sy) eq + ei = (2.27) 
where P is a polynomial in Sx and Sy of the form 

V = + .rf + ./erf) - f^S.^^.^ - kff^ 

-krflSx^ + ..., l,m = 1,2. (2.28) 

Resuming, we have derived the following system of equations 

Sl + Sl = Au (Jo) (2.29a) 

Sr = ^{Sx,Sy,Io,Ii,x,y,T) (2.29b) 

where we set 

LP (x, y, T, Sx,Sy) = --^-^'^ {Sx, Sy, X, y, r) - -^"^^^r"^- (2-30) 

In the construction above, an important role is played by the parameter 
with the condition < a < 1. It is now natural to ask whether such 
materials exist in reality. To provide with positive answer let us consider a 
medium which satisfies the following dispersion law 
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The formula (|2.31|) is referred to as the Cole-Cole dispersion law and it 
has been found experimentally by the Cole and Cole in 1941 5 . It is a 
phenomenoloigical modification of the "classical" Debye law (obtained from 
the formula H2.31() in the limit case a = 1). If eo ^-nd e depend on the 
intensity /, the function £i in the expansion (|2.19|1 is 

ei(/o,/i) = 7r(/o) h+e{Io)- (2.32) 

Thus, our general model is realizable in the Cole-Cole media. 

For sake of simplicity, let us consider a linearly polarized electric field 
E = e, where e is a constant unit vector and evaluate explicitly the 
nonlocal term D^r = e in the one-dimensional case 

/ + 00 
R{x' -x)e{l{x')) E{x') dx' (2.33) 
-oo 

If the nonlocal distribution function R {x' — x) is narrow one as it happens 
in the high frequency limit, one can expand e {I{x')) and E{x') around the 
point X. So, one gets 

Dn =Roe{I{x)) E{x)+Ri {e' {I{x)) E{x) + e{I{x)) E'{x)) + 



R2 { \e" {I{x)) E{x) + e' {I{x)) E'{x) + \e {I{x)) E"{x] 



where e' = dtjdx and e" = d'^e/dx'^ 

/+00 
R{x' -x) {x' -xY' dx' 
-oo 

is referred to as n— moment of the function R. If 



(2.34) 



Rn = ^, (2.35) 



in high frequency limit, one gets 

I?7V « ^ (r-o + inS, - \r2Sl + . . . ) e-^ a; ^ 00. (2.36) 

We note that narrower is the nonlocal distribution function R smaller are 
the higher moments. In this particular case the coefficients r„ are related 
one to another and they are expressed in terms of the fundamental quantity 
R. Thus, the equation of a;^~" order is 

5, = _|(,„ + i,,5._l.,S^ + ,..)_£!. (2.37) 
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An example of distribution R whose moments are of the form (|2.35|) is 
provided by the Gaussian distribution 

R{x) = ^ e-'^'''\ (2.38) 



In order to give a complete description of the physical system we should 
take into account the conservation law of the Poynting vector. The Poynting 
vector in complex representation is defined as follows [2]. 

P = EAB*. (2.39) 

In the more general case // 7^ 1, B should be replaced by H in the for- 
mula (E31I). 

The conservation law is (see pag. 35 in the reference [1]) 

V • P = 0. (2.40) 
Equation ()l.lb|) for time-oscillating fields looks like as follows 

VAE + icjB = 0. (2.41) 

Using the representation 

E = Eo e''^'^, B = Bo e'"^^ , 

one gets 

V AEo + iwVS'AEo = icjBo. (2.42) 
The leading order in high frequency limit gives 

Bo = V5AEo. (2.43) 

Thus, we have the following approximation of the Poynting vector 

P = Eo A Bo = Eo A (yS A ES 
(Eo • ES) VS - (Eo • V5) E5 

In virtue of the transversality condition ()2.22|) one finally has 

P = IVS. (2.44) 

where / = |Eo| . Note that P is parallel to the gradient of the phase VS. 
Using the expression 1)2. 44() in ()2.4U|) . one obtains 

VI -VS + IV^S = 0. (2.45) 
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Due to the slow dependence on the variable z, for phases of the form (|2.23j) 
and the asymptotic expansion / = Iq + uj^°'Ii + . . . , equation ()2.45|) becomes 

-^0 {Sxx + Syy) + IqxSx + loySy = (2.46a) 

h {SxX + Syy) + hxSx + hySy + kl^r = 0- (2.46b) 

The system (|2.46|) has to be considered together with the system for the 
phase ^TE^. 

The properties of present model have to be investigated by the joint 
analysis of the system ()2.29|) and the system ()2.46|) . In our discussions we 
will separate the compatibility analysis of the transverse equation (|2.29aj) 
and (|2.46a|) from the "evolution" equations (|2.29b|) and (|2.46b|) . 

3 Paraxial light beams. 
3.1 The nonlocal NLS equation. 

Many phenomenological models are based on the so-called nonlocal NLS 
equation |Hlini[Tni[niIIllI3IIllinilI7|. It describes paraxial light beams 
in nonlinear (and also nonlocal) media. In the present section we discuss 
the connection between the high frequency model discussed above and the 
nonlocal NLS equation. 

Let us consider a displacement vector of the form 

D = cjE + cj^D(3\ (3.1) 

where a is a constant parameter. The paraxial approximation in equa- 
tion ()2.3() is usually performed by the consideration of the slow variations of 
a small-amplitude electric field (or formally by the substitutions) 

(5,, dy, d,) ^ {adx, ady,a^d,) , E ^ aE e^^^ 

The leading nontrivial order in equation (|2.3() is 

2ia;— + ViE + L<j2D(3) = 0, (3.2) 
oz 

where V_l = {dx, dy). For a general nonlocal nonlinear medium we can write 
D(3) = f R{r'-r;a) N (/(r')) E(r') cfr' (3.3) 

where r = {x,y,z). We use the notation D^^^ to recall that, in the case of 
a local Kerr medium, the relation 1)3. 3|) is reduced to the cubic nonlinearity 
which is associated with the standard NLS equation j^j. The distribution 
i?(r' — r; a) characterizes the nonlocal response around the point r and a is 
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the "width" parameter (in the following it will be assumed to be depending 
on the frequency uj). N{I) is a certain nonlinear response depending on the 
intensity of the electric field / = |E| . 

We note that due to the paraxial approximation the nonlocal response 
along z can be neglected. As illustrative example to explain this fact let us 
consider a Gaussian nonlocal response 

i?(r',r)=^e-|r'-l^ (3.4) 



(2vr) 



2 



where |rp = x'^ + y'^ + z'^ . Due to the paraxial approximation we take into 
account of slow dependence on the variable z by the substitution z t~^z 

The Dirac 5— function in the left hand side of (|3.5|) implies that the response 
along the direction z becomes local. Thus, with the use of the distribu- 
tion (|3.5)) in the general definition of D^^^ the 3D integral is reduced to a 
2— dimensional one 

/+00 r+oo 
/ R{x'- X, y' -y,a) N (/(x, y, z)) E(x, y, z) dx dy. 
-00 J —00 

(3.6) 

The model considered above, can be further simplified under suitable as- 
sumptions on the widths of the nonlocal distribution R and the nonlinear 
response N{I). For sake of simplicity we focus on the 1 -|- 1— dimensional 
case 

(9E S^E 



where 



2iuj— + ^ + a;^D(3) = 3.7 

OZ OX"^ 

/ + 00 
R{x - x'; a) N {l{x')) B{x')dx'. (3.8) 

Let us define the width 5R of the nonlocal distribution R{x — x'; a) as the 
minimum such that 

R{x- x'; a) ~ 0, Vx' ^[x-5R,x + 6R] . (3.9) 

Of course, 6R depends on the width parameter a. Similarly, we can intro- 
duce the widths 5E and 6N of the electric field and the nonlinear response 
respectively. Suppose they satisfy the following conditions 

SR ~ 6N, SR « 5E. (3.10) 

Moreover, we assume the nonlinear response N[I) to be of the form 

N {I{x)) = N {X) , (3.11) 
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where X = ^ x and 7 oc 1/6N. Expanding E(x') and N {I{x')) in Taylor 
series around x, one gets the following approximation of the formula H3.8() 



d(3) ~ R{x -x';a) N {l{x'))dx'^B{x), (3.12) 

where we kept into account that due to equation (|3.11j) higher orders of 
the expansion of N (I(x)) are not negligible. Note that the formula ()3.12|) . 
in the case of nonlocal Kerr-type medium, leads to the nonlocal nonlinear 
Schrodinger equation discussed in the paper [TT] . 
For instance, given a bell-shape electric field E = Eg exp [— x^/2cr^ 

linear response of the form N{I) = I" = |Eo|^" exp — (7^)^ /2o"^ 

7 = \/2a, satisfies the condition (|3.1H) . Nevertheless, it is easy to see that 
only the validity of relations (|3.1U|) is sufficient to obtain the model ()3.12|) . 
For instance, for E = Eo/cosh^(x) and N {!) = 1°, condition (|3.10j) is 
verified for a large enough. Finally, we note that it is straightforward to 
generalize the previous considerations to construct a more general (2 + 1)D 
model. 



a non- 
where 



3.2 High frequency regimes. 

Now, we discuss the above in high frequency regime. 
With the representation 

E = Eo e^'^'^(^'^'^), 
the NNLS equation takes the form 

2S,+ (v^Sy = No{Io{x,y,z)), (3.13) 

where Vj_ = {dx, dy) and Nq^Iq) is the high frequency limit of the intensity 
law 

hm iV(/) = iVo(/o) 

(jJ—>00 

Paraxial approximation of the Poynting vector conservation law (|2.4()|) at 
the leading order on a gives the following equation on the xy— plane 

V±I -V^S + IVlS = 0. (3.14) 

An analysis of the (1 + 1)D reduction of the couple of the equations (|3.13j) 
and 1)3. 14(1 suggests that possible stable light beams (in this specific regime) 
exists only in (2 + 1)D. Indeed, the following example shows that the bell- 
shape initial beam profiles are no longer preserved. 
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If S in equations (|3.1,3j) and (|3.14|) does not depend on the variable y, one 
has 

25, + 52 = iVo (/o) (3.15a) 
IqSxx + IqxSx = 0. (3.15b) 

Integrating equation ()3.15b|) . one gets 

Sx = ^ (3.16a) 
-'0 

S, = l (^iVo(/o) - ^) , (3.16b) 
where 70 is an arbitrary real constant. In virtue of Sxz = Szx one obtains 
Jfiz + h — Ml 



where A'^^ := General solution of equation ()3.17|) . calculated by the 

characteristics method is provided by the following implicit relation 

^ + ^) = » (3.18) 

270 h J 

where $ is an arbitrary function of its argument. It is assigned by the initial 
profile of the intensity 

/o(x,0) = $-i(x). (3.19) 

It is well known ^H] that equations of type ()3.17|1 exhibit breaking wave 
phenomena for finite z. For example, a smooth initial profile, such that 
/q for X — > ±00 breaks at finite z where Iqx 00. 

Let us now consider a (2 + 1)D model for the following class of solutions 

S = kz + S{x,y,T), (3.20) 

where A: is a real constant and t = oj"z (with a > 0). A class of nonlocal 
perturbations with respect to the slow variable r will be discussed in the 
following subsection. Under these assumptions equations (|3.13j) and (|3.14j) 
assume the following form 

Sl + Sl = 4u{Io{x,y)) (3.21a) 
-^0 {Sxx + Syy) + IqxSx + loySy = (3.21b) 

where 4u(/o) = Nq{Iq) — 2k. We kept into account that due to the rule 

— =^"" — 
dz dr 
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2;— derivatives do not contribute in the limit to oo. We note also that Iq 
is the leading order term of the w""— asymptotic expansion of the intensity 

y, z) = Io{x, y, t) + LO~°'Ii{x, y,T) + ... . 

One refers to the functional dependence among u and Iq as intensity law. It 
is determined by the specific physical properties of the medium. Note that 
once u{Io) is given, the system 1)3. 21() is over determined. Its consistency 
condition will be discussed in the next section. 



3.3 Nonlocal perturbations. 

Let us consider a model for which the nonlocal nonlinear contribution to the 
displacement vector is of the following form 

D(3)=eE + D7v (3.22) 

where 

Djv = / R{r' -r,uj) N {l{r',z))E{r,z)dr' 

In particular, we are interested in the construction of a special class of weak 
nonlocal perturbations given by a nonlocal distribution function of the form 

M 

R{r-r',u)= ^ Ri,^ (r, .;) J^^'™) (r - r') , (3.23) 

l,m=0 

where r = (x, y) the distributions ij^''™) are defined by (|2.13|) and M is an 
arbitrary integer referred to as nonlocality degree. As a consequence one gets 

1=1 Lm=l 



Using the following asymptotic expansions on w " 

Eo(2;, y, z) = Eo(x, y, r) + a;""Ei(x, y, r) + 
s{x, y, z) = £0(2;, y, r) + w~"ei(x, y, r) + . . 



(3.25) 



and proceeding similarly to the previous section for a Cole-Cole medium ()2.31() . 
one gets from the NNLS equation in the leading orders oj"^ and uj^~°' the fol- 
lowing system 

^2 + ^2^ 4^^ (3.26a) 
Sr = ip{Sx,Sy, X, y, t) (3.26b) 

where Au = Eq — 2k and the function 99 is a polynomial in Sx and Sy. 
Equation H3.26a|) is the well known eikonal equation in two-dimensions. The 
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function (/? in equation (|.S.26b|) , in virtue of expression (|,3.24|) , is an M-degree 
polynomial in Sx and Sy and it describes an M-degree nonlocal response. 
Moreover, if we require that equation (|.S.26b|) possesses the phase inversion 
symmetry (5 —S) just like the eikonal equation, the function ip must 
contain only odd degree terms in Sx and Sy. We note also that the sys- 
tem of equations (|l-i.2(ij) has been derived first directly from the Maxwell 
equations [TB] , 

In local case we have 

ip = ip{z), u = u{x,y), (3.27) 

where if and u are certain function of their arguments. For a first degree 
nonlocality we have 

ip = aiSx + a2Sy, (3.28) 
and u must satisfy the following linear equation 

Uz = (aiu)^ + {a2u)y , (3.29) 

where ai and 02 are harmonic conjugate functions, that is they satisfy the 
Cauchy-Riemann conditions aix = a2y ctiy = —a^x- 



4 Transverse equations. 

The compatibility of the system of equations for the phase (|2.29|) and in- 
tensity 1)2. 46(1 results to be highly nontrivial problem. More specifically, one 
can see that it is not consistent for arbitrary intensity law. The construc- 
tion presented above suggest to separate the compatibility analysis on the 
xy— plane from the evolution along the z— direction. We will study the sys- 
tem of equations ()2.29a|) and (|2.46a|) for a general form of the intensity law 
and then we will discuss the "evolution" of the solutions with respect to the 
variable r. 

4.1 Elliptic intensity laws. 

Let us assume the function n(/o) to be invertible. Note that for various 
physically meaningful models such as the Kerr-type media 

u = ll, 7 > (4.1) 

and logarithmic saturable media 

u = log (1 + h/ht) (4.2) 

where the constant Iqi is the so-called threshold intensity, intensity law is 
the invertible one. 
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Then, we can consider the function /q = Io{u). We refer to it as inverse 
intensity law. Since 

V_l/o = I'o'V^u 

where Jq = dio/du, the system of equations (|2.29aj) and (|2.46aj) (or (|3.2H) ) 
looks hke as follows 

Sl + Sl = 4u (4.3a) 

-^0 {SxX + Syy) + I'q {UxSx + UySy) = 0. (4.3b) 

Substituting the expression for u from (|4.3a|) into (|4.3b|l . we obtain the 
following second order partial differential equation 

ASxX + BSyy + 2CSxy = 0, (4.4) 

where by definition Jq = log Jq and 

A = Jq'SI + 2, B = Jq'SI + 2, C = Jo'SxSy. (4.5) 

The equations of the form (|4.4|) are well studied in mathematics (see e.g. |19) 
and references therein). Their properties are critically depending on the 
signature of the discriminant 

A = AB- C^. (4.6) 

In particular, one can distinguish three cases, elliptic (A > 0), parabolic 
(A = 0) and hyperbolic (A < 0). In what follows, we will focus on the 
elliptic case motivated by the fact that the intensity law for Kerr-type 1)4.1(1 
and logarithmic saturable media ()4.2j) satisfy the ellipticity condition A > 
uniformly. We remark that these models are very relevant in practical 
applications since they are associated with the propagation of stable spatial 
solitons 1201 EH 

We would like to emphasize that elliptic second order nonlinear equa- 
tions of the form (|4.4|) possess several remarkable analytical and geometrical 
properties (see e.g. ^1 HSl IMl ES] ) • An interesting class of solutions of 
equation ()4.4|1 is provided by the Beltrami equation which is well known and 
well studied in the theory of elliptic systems of PDEs and in the theory 
of quasiconformal m,appinas\24:\ I25j . Let us introduce the complex variable 
z = X + iy {it should not be confused with real variable associated with 
the z— axis) and the complex gradient w = Sx — iSy. In these notations 
equation (|4.4() takes the form 

awz + bwz + aiBz + bwz = 0, (4-7) 

where 

a=^{A-B + 2iC) , b = ]^{A + B- 2iC) . (4.8) 
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It is immediate to verify that if w satifies the, so-called, nonlinear Beltrami 
equation 

Wz = tJ-{w,iI))wz, ij, = --. (4.9) 



then it is also a solution of equation (|4.7|) . It is easy to verify that ellipticity 
condition A > implies |^| < 1. Function fi possesses a remarkable geomet- 
rical meaning being the complex dilatation of the quasiconformal mapping 
w. 

Thus, the evolution of a light beam profile along the direction r turns 
out to be described by suitable deformations of quasiconformal mappings. 

More specifically, let us consider an input light beam profile such that 
/q 7^ inside a simply connected domain G and Iq = onto the smooth 
boundary F of G and outside it. Set uq = u{Iq = 0). Writing down the 
eikonal equation (|l'{.21a|) in terms of w and z, one gets 

w{z)id{z) = 4u. (4.10) 

From equation (|4.1U|) one gets {wl = for z £ T, that is F is mapped 

on the 2-yuo-radius circle. Assuming tt; to be assigned in such a way that, 
for instance, w{0) = 0, where S G, and the variation of the argument of 
the complex number w around F is Ap arg w = 27r, it can be proven that w 
is a homeomorphic mapping of domain G onto 2y^— radius disk F ' j22j . 
Of course, mapping w preserves the topology of domain G. In virtue of the 
assumption (|,'-{.2()|) the wavefront evolves according to the equation 

a;~°A;r + 5(x,y,r) = const. (4.11) 

We observe that the mapping w — Sx — iSy — {Sx, — Sy) can be also 
regarded as a two-dimensional vector field on the z— plane associated with 
transverse components of the wavefront normal unit vector 

n=^=(^ ^ UU) 

\VS\ VVA;2 + 4u' Vk^ + 4.u Vk'^ + ^uJ 

We remind that the vector n is, by definition, parallel to the Poynting vec- 
tor. This means that the mapping w encodes information about light-rays 
distribution around direction r. For example, let us consider the map- 
ping given in figure 14.11 Curves F and F' are oriented leaving domain on 
the right hand side. Under the assumptions mentioned above, there ex- 
ists a homeomorphism w of domain G onto G' acting in such a way that 
w{si,0) = (—1,0), w{0,S2) = (0,1) and it;(— S3,0) = (1,0). Consequently, 
the normal unit vector to wave-front on the boundary F is oriented in such 
a way that the light rays lie inside the rectangle circumscribing the domain 
G (see figure ITT|) . Recall that y-component of w reverse y-component of n. 
Then, this mapping describes a light beam 'trapped' around the direction r. 



17 



Conversely, if w(si,0) = (1,0), w{0,S2) = (0,-1) and w{-S3,0) = (0,-1) 
radiation spreads transversely far from r— direction. Note also that in 
both cases homeomorphism w is sense-reversing. If we consider a sense- 
preserving mapping such as, for instance, ?/;(si,0) = (1,0), w{0,S2) = (0,1) 
and w{—ss,0) = (—1,0), the beam spreads along x— direction and tends to 
be trapped along y— axis (see figure ESI- The boundary conditions and the 
orientation of the mapping allows us to control strongly the properties of 
the light beam. 

We emphasize, finally, that all of these observations can be generalized 
to the case of arbitrary n— connected domains. They could be used for 
description of interacting light beams. 




Figure 4.1: Sense-reversing mappings (solid line) describe beams which tend to 
be confined inside the rectangle or scattered (dashed line) outside it. 



Another class of solutions of equation 1)4. 7() can be obtained solving the 
equation 



bwz + awz = 0. 

Introducing the reciprocal coordinates by inversion of the system 

z = z{w, w) 
z = z(w, w). 



(4.13) 



(4.14) 




Figure 4.2: For the sense-preserving mapping the light beam spreads along 
X— direction and tends to be trapped along y— direction. 
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z w 




Figure 4.3: Homeomorphic mapping from the exterior of I?e on Dj^ such that T' 
goes onto F. Light rays associated with the points outside Dj^ are paraUel to the 
axis as much as e is smalL 

one converts equation ()4.13|) into the form of the linear Beltrami equation 

Zu, = iy{w,w)zw, (4.15) 

where i^{w, w) = —a/b and, due to ehipticity, < 1. The advantage of con- 
sideration of equation H4.13() in reciprocal coordinates is that the Beltrami 
equation is linear and it can be solved explicitly for different choices of the 
intensity law. Moreover, it is well known that several theorems in analytic 
functions theory can be rigorously generalized to the quasi-analytic functions 
which are the solutions of equation (|4.15)) (see e.g. Ref. In particular, 

a generalization of so-called Liouville theorem (Vekua's theorem) holds. If 
z = z{w, €)) is bounded on whole t/;— plane and satisfies the linear Beltrami 
equation ()4.15|) it can be shown that z{w, w) = constant '26'. We stress that 
for the constant solution, mapping from to— z— plane is singular and re- 
ciprocal transformation ()4.14l) is not defined. Then, any non-trivial solution 
of equation (|4.15|) must be singular somewhere on the complex plane and 
different type of singularities can occur, such as poles, essential singulari- 
ties, singularities of the derivatives etc. As illustrative example, we focus our 
attention on a solution z{w^ w) which have a simple pole at w = 0. As dis- 
cussed above, one can always consider a homeomorphism from D' = C \ D^^ 
to Dji, where is a disk of arbitrarily small radius e on the w— plane and 
Dfi is the i?-radius disk on the z— plane, mapping the boundary of on 
the boundary of D^. Inverse mapping w : — > D' , constructed in such 
a way, can be used to describe a beam "confined" around r— axis. Indeed, 
the transverse component of the vector n outside -D/j are arbitrarily small. 
Hence, the light rays can be settled down parallel to the r— axis with arbi- 
trary accuracy. So, the beam results to be self-guided around r as much as 
disk Df: is small. 

Coming back to nonlinear Beltrami equation (|4.9|) . we expect that for 
"mild enough" complex dilatations /i, Vekua's theorem still holds. In these 
cases, the only one bounded solution on whole 2— plane is w = constant. 
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In many physical situations, the intensity distribution on the z— plane, at 
certain r, goes to zero for z — > oo, or equivalently, one can says that inten- 
sity vanishes outside a big enough ii— radius disk Dr. Thus, outside Dr, 
refractive index assumes a constant value u = uq and the solutions of eikonal 
equation (|3.21a|l is 

5 = cqx + ciy + C3, (4.16) 

where cq, ci and C3 are constants and the condition c^ + c\ = 4uo holds. For 
a paraxial beam we have cq = ci = 0. As the consequence w — Sx — 
Co — ici = and, of course it satisfies Beltrami equation in C \ Dr. In 
virtue of Vekua's theorem, the only one bounded solution is u; = 0. Then, 
any non-trivial solution must be singular somewhere on the plane. In our 
example, wavefront is approximately plane for z E C \ Dr and possesses 
singularity inside Dr. 

In the next subsection we will discuss the so-called optical vortices ob- 
tained from equation ()4.4() assuming the phase to be a harmonic function on 
the xy— plane. We will see that, in this case, equation ()4.4() is equivalent to 
the minimal surfaces equation. We refer to this class of solutions as minimal 
sector in order to distinguish it from the Beltrami sector which is associated 
with the solutions of equations (|4.9j) and (|4.15j) . 



4.2 Helicoidal wavefront s. 

The successive approximations method is a general approach to solve the 
linear Beltrami equation (it has been demonstrated by Tricomi |27j ) . Never- 
theless, calculation of exact 'explicit' solutions can be the difficult task and 
the chances of success strongly depend on the form of the complex dilata- 
tion. Incidentally, we note that solutions possessing cylindrical symmetry 
such that 

S = S{r), u = u{r), r = ^/x"^ + y"^, (4-17) 

are not compatible with intensity law. It is straightforward to verify that 
equations (|3.2H) along with assumptions (|4.17l) imply that intensity depends 
explicitly on z— axis distance r = \/x2 + y'^ 

where c is an arbitrary constant. 

Here, we will consider the solutions of equation (|4.4j) connected with 
the, so-called, minimal surfaces. The minimal surface equation looks like as 
follows (see e.g. Ref. [2H1) 

(1 + Si) SxX + (1 + Si) Syy - 2SxSySxy = 0. (4.19) 

We restrict ourselves to the class of solutions of equation (|4.4|) which are 
also harmonic on the xy— plane, that is 

SxX + Syy = 0. (4.20) 
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Using the condition (|4.2fl|) in equation (|4.4j) , one gets the equation 



S^Sxx + ^ySyy ~\~ '^SxSySxy — 0, (4.21) 

for any elhptic intensity law. It is straightforward to check that equa- 



tions (|4.21j) and (|4.19j) coincide for harmonic solutions. In other words, 
a class of solutions of equation 1)4. 4|) . whatever is the inverse intensity law 
-^0 = loiu), is just given by the class of the harmonic minimal surfaces. An 
important result for us is that the only non-trivial harmonic minimal surface 
in Cartesian coordinates is the helicoid ^1. It can be written as follows 



5 = Karctan( ) , (4.22) 



where K is an arbitrary constant. One can check that the function S given 
by (|^^ satisfies equations ^J^, and (|^^ simultaneously. Equa- 

tion of corresponding wavefronts is 

S = uj'^T + K arctan f = const, (4.23) 



where the constant K is the "pitch" of the helicoid. In particular, the 
expression 1)4. 23() describes the edge-screw dislocations discussed first time 
experimentally by Brynghdal in 1973 29^ and theoretically by Nye and 
Berry jHUI in 1974. We highlight that helicoidal wavefronts exist in both 
linear and (nonlocal) nonlinear regimes. They are associated with singular- 
ity of the phase (phase defects) which appears as topological defects of the 
interferograms. These class of phase defects has important phenomenologi- 
cal consequences connected to the optical vortices EH ESI EH (see also 
Refs. 133 EEl and references therein). If we assume, for simplicity, (3 = 
in 1)4. 22() . one has the pure screw dislocations. The complex gradient asso- 
ciated with the helicoid 

K 

w = -i— (4.24) 

z 

is a meromorphic function on z with the simple pole singularity at the origin. 
Nevertheless, normal vector to the wavefront, whose components coincides 
(up to a sign) with real and imaginary parts of w, is not defined. Indeed, 
vector w = {Sx, —Sy) has no limit as {x,y) — > 0. Figure ^31 shows examples 
of helicoidal wavefronts 1)4. 23() usually parametrized as follows 

x = vcost, y = vsmt, z = Kt. (4-25) 

One-start helicoid shown in figure ITlk is obtained for < v < +oo; two- 
start helicoids shown in figure I4.4b .c are obtained for — cxd < t; < -|-cxd. 
Refractive index 

u = ^ (4.26) 
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Figure 4.4: Helicoidal structures of wavefront around r— axis: a) one-start 
right-screw {K = 1); b) two-start right screw {K = 1); c) two-start left 
screw {K = —1). 



has cylindrical symmetry around r— axis and displays divergence at z = 
0. This means that around the origin geometric optics approximation fails 
and wave effects become relevant. In particular, necessary condition for 
the existence of singular wavefronts is that intensity vanishes where phase 
function is singular 130J . Indeed, in this region the interference phenomenon 
is no more negligible and can realize this condition. 

It is interesting to evaluate the effect of the present class of nonlocal 
perturbations 

Sr = f{Sx,Sy,X,y,T) 

discussed in the sections HI and 121 on the pure (/3 = 0) helicoidal wave- 
front (|4.22|) . In particular we note that the solution (|4.22j) of the equa- 
tion (|4.21|) is defined up to an additive arbitrary function of r— variable 
■0(t), such that ip' = '^{t). Then, one has 

S = K arctan (y/x) + ?/'(t). 

Let us observe that the function tpi^r) is specified by the function ip evaluated 
on the helicoid. For example, in the case of the first degree nonlocality one 
has (see equations p. 28(1 and (|3.29|) ') 

iP'{t) = aA + a2Sy 
Ur = {aiu)^ + {a2u) . 

Evaluating equations (|4.27|) on the helicoid one gets 

(x^ + y^) aix = xai ' 
-yai + xa2 = fir) (x^ + y^) . 



2 , 2^ (4-28) 



First degree nonlocal deformations of the helicoid exist only for the nonlocal 
data ai and 02 satisfying the system ()4.28|) . For example, trivial solutions 
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ai = a2 = (p = corresponds to the local case discussed above. A simple 
non-trivial solution ai = —^y, 02 = jx and ^' = 7, where 7 is an arbitrary 
constant, provides us with the following wavefront 

S = uj'^T + ^ — arctan (-) = const. (4.29) 

For 7 = equation 1)4. 29(1 coincides with equation 1)4. 23(1 . As a consequence 
of nonlocal response the helicoid's pitch is compressed if 7 > and stretched 
if 7 < 0. 



5 Nonlocal perturbations and dVN hierarchy. 

In this section we will discuss a particular class of nonlocal perturbations 
which are connected with an integrable hierarchy of PDEs, namely, the 
dispersionless Veselov-Novikov hierarchy. The methods to solve these equa- 
tions, in particular, the reduction method based on the symmetry constraints 
and the quasiclassical 5— dressing approach will be considered in sections 6 
and 7. 

Using the complex variables 2; = x + iy, z = x — iy, one rewrites equa- 
tions (|2.29|) as follows 

SzSz = u{z,z,t) (5.1) 
Sr = p{Sz,Ss,z,z,t) . (5.2) 

The compatibility condition of equations ()2.29|) imposes constraints on the 
possible forms of the function ip, namely 

SzPz + SziPz + Uzif' + Uzip" = Ur, (5.3) 

where 

oSz oSz 
The simplest cases have been already discussed above and corresponds to 
the local and first degree nonlocal responses, given by the function ip of the 
forms (|3.27|) and 1)3. 28() respectively. The quadratic case 

if = aSl + aSl + l3Sz + PS, + 7 + 7, (5-5) 

is equivalent to the linear one. Indeed, the consistency of equations (|5.3jl 
and (|5.1|) implies 

a = 0, /3 = (3{z), 7 = const, 

and the (|^ is reduced to the (jSUHl)- 
The cubic function 

if = aSl + aSl + l3Sl + j3Sl + 7^^ + ^S-z + 5 + ~5 (5.6) 



23 



obeys equations (|5..S|) and (|5.1j) if 

a = a{z), /5 = 0, 
Iz = -azU - 3auz, 

6 = const, (5-7) 

and one has the equation 

Ur = (ju)^ + {ju)^ . (5.8) 

In the particular case a = 1 and, consequently 72 = —3uz, equation (|5.8j) is 
the dispersionless Veselov-Novikov (dVN) equation introduced in P71 

In a similar way, repeating the procedures for higher degree polynomials, 
one can construct higher order nonlocal perturbations. 

Besides, if we formally admit all possible degrees of Sz and Sz in the right 
hand side of equation (|5.2() . one has an infinite family of nonlinear equations, 
governing the nonlocal deformations of the wavefronts and "refractive index" 
u. 

Moreover, the polynomial dependence of these deformations should be 
compatible with certain constraints. The request that equation (|5.2j) re- 
spects the symmetry for phase inversion, S — > —S, of the eikonal equa- 
tion gives the hierarchy of nonlocal deformations of the form 

n n 

J2 u^Sl^-^ + J2 ^mS?'"-^ (5.9) 

m=l m=l 

Note that the constant terms which have appeared in (|5.8|) is. in fact, ir- 
relevant. Hence, these polynomial deformations assume the general form 

n 

Sr=Yl i^mS'^z"''' + , n = 1, 2, 3, . . . (5.10) 

m=l 

where Um are certain functions on u. 

In the = 1 one gets the dVN equation mentioned above (n = 2) 

and the so-called dVN hierarchy of nonlinear PDEs. This means that the 
dVN hierarchy is associated with a specific class of "integrable" nonlocal 
perturbations. In particular, there is a one-to-one correspondence between 
nonlocal degrees and the equations of the dVN hierarchy. It is straightfor- 
ward to verify that the helicoidal wavefronts are not preserved under dVN 
hierarchy nonlocal perturbations. 

6 Integrable systems. 

The celebrated Korteweg-de-Vries (KdV) equation 

Ut + 6uU:j: + Uxxx = (6.1) 
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has been introduced to describe long, one-dimensional, small amplitude shal- 
low water wave propagation j39j . By a numerical analysis Zabusky and 
Kruskal [lOj observed main feature of the solitons, that is they "pass one 
another without losing their identity" . Solitons, can be seen as a purely 
nonlinear phenomenon, where the linear dispersion is "balanced" by the 
nonlinearity. Fundamental step has been taken in 1967 with the famous 
paper by Gardner, Green, Kruskal and Miura, with the introduction of the 
so-called inverse spectral transform method (1ST). The 1ST approach works 
as a nonlinear Fourier transform according to the following scheme 

u{x,0) Scattering data at t = ^ evolution 
Scattering data at t = — > u{x,t'). 

By means of the spectral transform the initial profile u{x, 0) is associated 
with a set of spectral data at the time t = 0. In the space of spectral data the 
system results to be linearized and the evolution is trivial. Then, once the 
spectral data at the time t = t' are obtained, inverse spectral transform gives 
the evolved profile u{x,t'), solution of nonlinear PDE. The method is based 
on the representation of a nonlinear PDE as the compatibility condition of 
a pair of linear problems. In the KdV case the linear system is 

- - u{x,t)ip = X^il;, (6.2a) 

where A is usually called spectral parameter. The system 1)6. 2() can be re- 
garded also as the compatibility condition of the couple of operators {Lax 
pair) 

[^1,^2] =0 (6.3) 



where 



52 

Li = ^ +n(x,t) + (6.4) 

The nonlinear Schrodinger (NLS) equation 

iqt + Qxx + "^s Iql"^ q = (6.6) 
and sine-Gordon (SG) equation 

utt - Uxx + smu = (6.7) 
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are other remarkable examples of (1 + 1)D nonlinear integrable PDEs. NLS 
equation has applications, for instance, in nonlinear optics jl] and Bose- 
Einstein condensation [2] and SG equation appears in geometry to de- 
scribe negative curvature surfaces j42j, in physics in the study of flux prop- 
agation in Josephson junctions, nonlinear optics and quantum field the- 
ory (see e.g. [Sni)- Thereafter, the method has been generalized to (2 -|- 
1)— dimensional soliton equations, such as the Kadomtsev-Petviashvili equa- 
tion 

3 3 

2 4 ^ (6.8) 

Vx = Uy. 

Linear problems for (|6.8j) are 

ipy = ipxx + uip (6.9a) 

3 /3 3 \ 

ipt = ipxxx + -^uipx + ( -jUx -f j V = (6.9b) 

where Vx = Uy. Equation ()6.8() is the {2 + l)D generalization of KdV describ- 
ing weakly two-dimensional long, shallow water waves. Another (2 -|- 1)D 
integrable generalization of KdV equation is given by the Nizhnik-Veselov- 
Novikov (NVN) equation 

Ut + CiU^^^+C2Urfr^ri-'iCi(^d'^^Ur^ - 3C2 {ud~^U^) ^ = (6.10) 

where = x + sy, rj = x — sy. Equation (|6.10|) has been introduced by 
Nizhnik in the case s = 1 jHj and by Veselov and Novikov in the case s = i, 
ci = C2 = 1 ^ . Unlike the KP equation no physical applications of the NVN 
equation is known at the moment. The Veselov-Novikov (VN) equation 
{z = X + iy) 

ut = {uV)^ + {uV). + Uzzz + Uzzz (6.11a) 
V-z = -2,uz, (6.11b) 

is equivalent to the compatibility condition for equations 

tpt = Ipzzz + i'zzz + (Vlpz) + (Vlpz) . 

Note also that the local (2 -|- 1)D NLS equation is not integrable by the 1ST 
method. 

Recently the, so-called, quasi-classical limit (or dispersionless) limit of 
soliton equations has attracted a great interest (see e.g. [ini)-[SH]- Disper- 
sionless limit of soliton equations is performed formally introducing a slow 
variable expansion by the substitution 

tn ^ tn/e (6.12) 
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and introducing the asymptotic expansion of the form 

u(^^^ =u{tn) + 0{e), e^O. (6.13) 

where tn the set of the "time" variables. For example, in the KdV equation 
case one identify ti = x, t2 = y, while for KP equation it is h = x, t2 = y 
and ts = t. 

In terms of the linear problems, this procedure corresponds to the fol- 
lowing representation of the function ■0 

s 

where 

^(a;^) =5(t„) + 0(e). 

For example, in the leading order equation (|6.8() is the dispersionless KP 
(dKP) equation 



3 3 

Ut = -UUx + -Vy 

2 4 ^ (6.14) 



Vx = 



and the linear problems 1)6.9(1 becomes the Hamilton-Jacobi equations 



Sy ^ X ^ 

a o3 3 ^ 3 3 (6-15) 
St = + -uSx + -Ux + -V. 

The dKP equation is known in physics as Zabolotskaya-Khokhlov equa- 
tion |59] . 

In the following, in connection with the model of nonlocal nonlinear 
optics discussed above, we will focus on the dispersionless Veselov-Novikov 
(dVN) equation. It looks like as follows 

Ut = {uV)^ + iuV)- 
Vg = -3uz. 

It is obtained as the compatibility condition of the following system 

, , - (6.17) 

St = Si + Si + VS, + VS-,. 

Dispersionless systems are of interest for both physical and mathematical 
reasons. They have applications in hydro dynamincs, magneto- hydrodynamics, 
Laplacian growth [201. They arise also in the framework of the Whitham 
averaging method for calculation of small amplitude modulations of soliton 
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equations solutions |18[ 160] . Moreover, we will see in the following that the 
dVN equation and the hierarchy associated are relevant in the description 
of specific high frequency regimes in nonlocal nonlinear optics. 

Different approaches can be used to study dispersionless systems. Here 
we will consider the reduction method based on the symmetry constraints 
and the 5— dressing method. In particular, their applications to the eikonal 
equation and the nonlocal perturbations associated with the dVN equation 
will be also discussed. 



7 Symmetry constraints. 

Here we will recall shortly the definition of symmetry constraint and discuss 
its specific application to the dVN equation as an approach for construc- 
tion of (1 + ID) hydrodynamic type reductions. Let us consider a partial 
differential equation for the scalar function u = u{t) = u {ti,t2, ■ ■ ■) 

F {u,ut^,ut,t,,---) = 0, (7.1) 

where = du/dti. By definition, a symmetry of equation ()7.1() is a trans- 
formation u{t) u'{t'), such that u'{t') is again a solution of (|7.H) (for more 
details see e.g., [HH). Infinitesimal continuous symmetry transformations 

t'i = ti + 6ti; u' = u + 6u = u + eue (7.2) 

are defined by the linearized equation 1)7. 1|) 

L5u = 0, (7.3) 

where L is the Gateaux derivative of F 



dF ( d , , , 

Lou := ( ^ + ^'"o ^ + ^^e) ,•••)>■ 



(7.4) 



e=0 



Any linear superposition 6u = CkS^u of infinitesimal symmetries 6kU is, 
obviously, an infinitesimal symmetry. By definition, a symmetry constraint 
is a requirement that certain superposition of infinitesimal symmetries van- 
ishes, i.e. 

'^CkSkU = 0. (7.5) 

k 

Since null function is a symmetry of equation (|7.H) . the constraint (|7.5)) is 
compatible with equation ()7.1() . Symmetry constraints allow us to select a 
class of solutions which possess some invariance properties. For instance, 
well-known symmetry constraint 6u = eut^, = 0, selects solutions which are 
stationary with respect to the "time" tk- 



28 



7.1 The dKP equation. 

The linearized equation associate with the dKP equation l|6.14|) is 

3 3 

{6lo), = {6u)y. (7.6) 
Solutions are infinitesimal symmetries of dKP. 

Theorem 1 Suppose Si and Si are arbitrary solutions of the Hamilton- 
Jacobi equations ()6.17() . Then the quantity 

N 

5u = '^Ci (^Si - Si^ , 
1=1 

where Ci are arbitrary constants, is a symmetry of dKP equation. 



Proof. It is straightforward to check that i Si — Si) satisfies equa- 

V /XX 

tion dZSI). □ 

This type of symmetries has been introduced for the first time in j63j . 
within the quasiclassical 9-dressing approach. As a simple example, let us 
consider the following symmetry constraint 

Ux = Sxx- (7.7) 

Under this constraint the Hamilton-Jacobi system (|(i.l5|) gives rise to the fol- 
lowing hydrodynamic type system (the dispersionless nonlinear Schrodinger 
equation, see (62)): 

Uy = [u^ + u)x, 

uy = 2 {uu), (7.8) 

where u = dSx/dX. 

7.2 Real dVN equation 

Let us focus on the case of real-valued u. Infinitesimal symmetries 5u of the 
dVN equation obey the equations 

{5u)^ = {V6u + u6V)^ + {V6u + 6Vu). 

V, = -3u,; {6V), = -3{6u)^. (7.9) 

Theorem 2 Suppose Si and Si are solutions of the Hamilton-Jacobi equa- 
tions (|6.17|1 . Then the quantity 

N 

6u = Y,Ci[Si-Si)^_, (7.10) 
1=1 

where Ci are arbitrary constants, is a symmetry of the dVN equation. 
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Proof. It is straightforward to check that (^Si — Si^ _ satisfies equa- 
tion dZSI). □ 

In particular, one can choose Si = S{X = Aj) and Si = S{X = + fii). 
In the case /ij ^ and q = Ci/fii, one has the class of symmetries given by 

N 

5u = ^^Ci4>i-^z (7.11a) 

i=l 

0^ = ^(A = A,). (7.11b) 

In what follows we will discuss three particular cases of real reductions, 
providing real solutions of the dVN equation. 

If 5 is a solution of Hamilton- Jacobi equations (|6.17|1 . then —5 is a solu- 
tion as well. Thus, for real- valued S {S = S), specializing constraint H7.1U() 
for = 1, we have a simple constraint 

u, = iS),,. (7.12) 

Let us introduce the functions pi := Sx and p2 := Sy. Thus, the symmetry 
constraint ()7.12|) can be written as follows 

= ^ {Sxx + Syy) = ^ {pix + P2y) ■ (7-13) 

In order to analyze constraint 1)7. 13() it is more convenient to consider equa- 
tions (|6.17|) in Cartesian coordinates (x, y), i.e., 

Sl + Sl = Au (7.14a) 

St = \sl - \SxSl + ViSx + V2Sy, (7.14b) 
where V = Vi + iV2, while the dVN equation acquires the form 

ut = (uVi)^ + {uV2)y (7.15a) 

Vix - V2y = -?>Ux (7.15b) 

V2x + Viy = 2,uy. (7.15c) 

Substituting (|7.14a|) in (|7.13j) . one obtains the following hydrodynamic type 
system 

Pi \ _ ( I \ [ pi 



P2 J y \2pi-l 2p2 J \ P2 J ^ 

Now, let us focus on definition Vz := —2>Uz. Differentiating it with 
respect to x, using constraint ()7.12|) and equations 1)7. 16() . one obtains the 
equations 

Vix = -^Pix + ^{pI + pf)^ 
3 

y2x = -^P2x, 
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which can be trivially integrated providing the following explicit formulas 
for Vi and V2 in terms of pi and p2- 



V2 = Ip2. (7.17) 

At this point we can derive t-dependent equations for pi and p2- Differenti- 
ating equation ()7.14b|) and using (|7.16|) and l|7.17|) . one obtains the system 




(7.18) 



where 

^11 = 3/91 (pi - 1) , A12 = 3p2, 
A21 = 3p2 (2/91 - 1) , A22 = 3pi (pi - 1) + 6pl 

Common solutions {pi,p2) of the systems (|7.16() and 1)7. 18() provides us with 
the solution u = [pf + pf) /4 of the dVN equation (|7.15|) . 

7.3 Admissible intensity laws. 

The analysis of the compatibility between the class of nonlocal perturbations 
associated with the dVN hierarchy and the intensity conservation law (|3.21b|) 
shows that there are no nontrivial solutions for an arbitrary form of the 
intensity law. In particular, intensity law appears to be a quite restrictive 
constraint making the problem of compatible nonlocal responses rather non- 
trivial. Nevertheless, it is possible to see, by an explicit example, that there 
exists nontrivial intensity laws such that the system (|3.26j) and the intensity 
conservation law (|3.21b|) are compatible with the dVN hierarchy. Here, as 
illustrative example, we focus on the third degree of nonlocality associated 
with the dVN equation. 

In order to do that, we consider hydrodynamic type reductions of dVN 
equation. They have been found using symmetry constraint of the form 

V]_S = Ux 

or equivalently Ux = Szz- Under such a constraint one gets the hydrody- 
namic type system (|7.1()|) . (|7.18|) . 

Looking for solutions of the system 1)7. 16() and 1)7. 18() such that p2 = 
P2 {pi ) ) one finds that pi and p2 are given implicitly by the following algebraic 
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system 



x + G'y + H'z - $ {pi 



+ 



2c - log ( g + \/l + 



P2 = 2 



q=P2± \ P2 + 2pi - 1 



(7.19) 
(7.20) 



where c is an arbitrary constant, G = p2{pi), H = p\ — ^pf + ^P2{pi) and 
'prime' means the derivative with respect to pi. 

Differentiating eikonal equation (|3.26a|) with respect to x and taking into 
account that Ux = V^S", one gets 



V 



Sx 



v±s + vls = 0. 



(7.21) 



Comparing equation (|7.2H) with intensity conservation equation ()3.21b|) . 
which can be written equivalently as follows 



V^S-V^ (log/o)+Vi5 = 



(7.22) 



one gets the following simple relation among intensity and pi— component 
of the gradient 

/o = Ce-'^, (7.23) 



where C is an arbitrary real constant. Finally, eikonal equation provides us 
with intensity law 



u{Io) = log 



C 



+ 



1 



P2 -2 log 



C 



(7.24) 



where last term in r.h.s. is given by the algebraic relation 1)7.20(1 . 



8 The quasiclassical (9-dressing method. 
8.1 The eikonal equation. 

In this and next sections we will demonstrate that the plane eikonal equation 
and dVN hierarchy both for complex and real refractive indices are treatable 
by the quasiclassical 3-dressing method. 

The quasiclassical 5-dressing method is based on the nonlinear Beltrami 
equation f531 

S-^ = W{\~X;Sx) (8.1) 
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where ^(A, A) is a complex valued function, A is the complex variable S\ = 
dS/dX and W (the quasiclassical 3-data) is an analytic function of Sx 



W{X,\Sx)=J2'^n{X,X){Sxr, (8.2) 



n=0 

with some, in general, arbitrary functions U)„(A,A). 

To construct integrable equations one has to specify the domain G (in 
the complex plane C) of support for the function W {W = 0, X £ C/G) 
and look for solution of (|<S.1|) in the form S = Sq + S, where the func- 
tion Sq is analytic inside G, while S is analytic outside G. In order to 
construct the eikonal equation on the plane we choose G as the ring 
P = {a G C : i < |A| < a}, where a is an arbitrary real number (a > 1), 
and select B — type solutions satisfying the constraint 

S [-X,-X) = -S {X,X) . (8.3) 

Then, we choose 

So = zX + ^. (8.4) 
A 

Due to the analyticity of S outside the ring and the property (|8.3|) one has 

oo ct(oo) 

^ = E#Tr; ^-oo, (8.5) 

n=0 

and 

oo 

5 = 5^5£ViA'"+' A-0. (8.6) 

n=0 

In particular, 5 (0, 0) = 0. 

An important property of the nonlinear (9-problem (|8.H) is that the deriva- 
tives f = St oi S with respect to any independent variable t, obeys the linear 
Beltrami equation 

f-x = W'{X,X;Sx)fx (8.7) 

where W' (A, A; (p) = (A, A; (p) . Equations (|8.7j) has two basic properties, 
namely, 1) any differentiable function of solutions fi, ■ ■ ■ , fn is again a solu- 
tion; 2) under certain mild conditions on W' , a bounded solution / which is 
equal to zero at certain point Aq € C, vanishes identically (Vekua's theorem) 

121. 

These two properties allows us to construct an equation of the form 
^iSz,Sz) = with certain function Q. Indeed, taking into account (|8.4|) . 
one has 

Sz = X + Sz, (8.8) 
S-z = \ + S,, (8.9) 
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i.e. Sz has a pole at A = oo, while S'g has a pole at A = 0. The product SzSg 
is again a solution of the linear Beltrami equation 1)8. 7|) and it is bounded 
on the complex plane since 

SzS-z = l + ^Sz + XSz + SzS-z (8.10) 

where 5^ (A = 0) = and XSz s[^^ as A ^ oo. 

Subtracting l-\-S[^^ from the r.h.s of equation 1)8. 1U() . one gets a solution 
of equation (|8.7|) which is bounded in C and vanishes as A — > oo. According 
to the Vekua's theorem it is equal to zero for all A. Thus we get the equation 

SzS-z = u{z,z) (8.11) 

where 

U = 1 + S['^J. (8.12) 

In the Cartesian coordinates x, y defined hy z = x + iy, equation (|8.11|) is 
the standard two-dimensional eikonal equation discussed above 

{V±Sf = Au (8.13) 

where we recall that ~ 

Using the 5-problem (|8.H) . one can, in principle, construct solutions 
of equation (|8.11|) . So, the quasiclassical 5-dressing method allows us to 
treat the plane eikonal equation 1)8. 11(1 in a way similar to dKP and d2DTL 
equations. We note that the phase function S in (|8.1H) depends also on 
the complex variables A and A. Curves S{z, z) = const define wavefronts. 
The 9-dressing approach provides us also with the equation of light rays. 
Indeed, since r.h.s. of 1)8. 13(1 does not depend on A and A, the differentiation 
of (|8.13|1 with respect of A (or A) gives 

V5-V(/> = (8.14) 

where (j) = S\ (or (p = S-^)- So, the curves S = const and (p = const 
are reciprocally orthogonal and, hence, the latter ones are nothing but the 
trajectories of propagating light. Thus, the 9-dressing approach provides us 
with all characteristics of the propagating light on the plane. Note that any 
differentiable function (f>{S\,S^) is the solution of equation (|8.14|1 too. 

In general, within the 9-dressing approach one has a complex-valued 
phase function 5 and, consequently, a complex refractive index. To guaran- 
tee the reality of u, it is sufficient to impose the following constraint on S 

S{X,X)=s(^j,j^ . (8.15) 

Indeed, taking the complex conjugation of equation (|8.11|) . using the differ- 
ential consequences (with respect to z and z), of the above constraint and 
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taking into account the independence of the l.h.s. of equation (|8.11|) on A, A, 

one gets 

u (xn) = S, (A, A) (A, A) = S, = u (x„) , (8.16) 

i.e. the "refractive index" u is real one. The constraint (|8.15|) leads to the 
relations 5'2n+i ~ ^2n+i- Moreover, this constraint implies also that the 
function S is real-valued on the unit circle |A| = I {S (A, A) = 5 (A, A), 
|A| = 1). This provides us with the physical wavefronts. 

The (9-approach reveals also the connection between geometrical optics 
and the theory of the, so-called, quasiconformal mappings on the plane. 
Quasi-conformal mappings represents themselves a very natural and impor- 
tant extension of the well-known conformal mappings (see e.g. HHI)- In 
contrast to the conformal mappings the quasi-conformal mappings are given 
by non-analytic functions, in particular, by solutions of the Beltrami equa- 
tion. 

A solution of the nonlinear Beltrami equation (|8.H) defines a quasi- 
conformal mapping of the complex plane C |24[ I25j . In our case we have a 
mapping S (A, A) which is conformal outside the ring G and quasi-conformal 
inside G. Such mapping referred as the conformal mapping with quasi- 
conformal extension. So, the quasi-conformal mappings of this type which 
obey, in addition, the properties (|8.3() . ()8.4() and 1)8. 15(1 provide us with the 
solutions of the plane eikonal equation (|8.11|) . In particular, wavefronts given 
by S (A, A; z, z) = const are level sets of such quasi-conformal mappings. In 
more details, the interconnection between quasi-conformal mappings and 
geometrical optics on the plane will be discussed elsewhere. 

8.2 The dVN hierarchy. 

In this section we will apply the quasiclassical 9-dressing method to the dVN 
hierarchy. For this purpose we consider again the nonlinear 9-problem (|8.1|) , 
(|8.2|) on the ring G with the constraints (|8.3|) and (|8.15j) and introduce 
independent variables Xn, x^ via 

oo oo 

So = Y. x„A2«~i + ^ x„A-2"+i. (8.17) 

n=l n=l 

The derivatives Sx„ and Sx„ obey the linear Beltrami equation (|8.7|1 . and 
using the Vekua's theorem one can construct an infinite set of equations of 
the form 

Q{Xn,Xn,Sx„,SxJ = 0. (8.18) 

Repeating the procedure described in the previous section, one gets the 
plane eikonal equation for z = xi. In a similar manner, for the variable 
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T = X2 = X2, taking into account that Sr = + ^ + Sr, one obtains the 
equation 

Sr = Si + Si + VS, + VSs (8.19) 

where V = —3S[^'^ = —Sd^^Uz- Evaluating the terms of the order in 
the both sides of equation H8.19() , one gets the dVN equation 

Ur = -3 (ud^^u,)^ - 3 {ud;\,). . (8.20) 

In similar manner, setting r = X3 = X3 one has 

^ = Sl + Sl + asSi + asSi + a^S, + a^S-, (8.21) 

and consequently the equations 

Ur = {a-iu)^ + (ain)j (8.22a) 
+ asg = (8.22b) 
ai2 + ua^z + 3a3U^ = (8.22c) 

Considering the higher variables r„ = x„ = x„, (n = 2,3,...), one con- 
structs the whole dVN hierarchy. It is a straightforward to check that the 
constraint (|8.15j) is compatible with equations (|8.19|) and higher ones. So, 
the 0-dressing scheme under this constraint provides us with the real- valued 
solutions of the dVN hierarchy. 

If one relaxes the reality condition for u (i.e. does not impose the con- 
straint (|8.15j) ). then there is a wider family of integrable deformations of 
the plane eikonal equation with complex refractive index. To build these 
deformations one again considers the 9-problem (|8.1|) . chooses the domain 
G as the ring (as before), imposes the 5-constraint (|8.3|) . but now chooses 

as follows 

5o = ^ x„A2"-i + ynA"2"+i (8.23) 

n=l n=l 

where yi = xi = z. Repeating the previous construction, one gets again the 
eikonal equation (|8.11|) . but now with a complex- valued u. Considering the 
derivatives Sx^^ Sy^ with n = 2, 3, . . . , one obtains the two set of equations 

n 

Sxr, = Y^Um{Xn,yn){Szf''^~^-, n = 2, 3, . . . (8.24) 

m=l 
n 

Syn=^^Vm{xn,yn){Szf''"'^^ , n = 2, 3, . . . (8.25) 

m=l 

Equations H8.11|) and H8.24|) give rise to the hierarchy of equations 

Ux^= Fn{u,Uz,Uz) , (8.26) 
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the simplest of which is of the form 

u^^ = -3 [ud^^u,]^ (8.27) 

where 

S:,, = Sl + VS,. (8.28) 

Equations 1)8. and (|8.25p generate the hierarchy of equations Uy„ = F^, 
the simplest of which is given by 

Uy, = -?,{ud-^u,)., (8.29) 

and 

Sy2 = Si + VS-,. (8.30) 

For both of these hierarchies the quantity J udz A dz is the integral of 
motion as for the dVN hierarchy. It is easy to see that equations 1)8. 27() 
and (|8.29j) imply the dVN equation ()8.2n|) for the variable t = {x2 + 2/2) /2. 

8.3 Characterization of 5-data. 

As we have seen, the constraints (|8.3)) and ()8.15|) guarantee that one will 
get the eikonal equation ()8.11|) with real valued u. In this section we will 
discuss the characterization conditions for 9-data W (A, A, 5a) which provide 
us with such result. In general, if one considers 3-problem with a kernel 
defined in a ring and the function Sq singular at two points, e.g. A = 
and A = 00, without i?-constraints (|8.3|) . one constructs the dispersionless 
Laplace hierarchy [HHl associated with the quasiclassical limit of the Laplace 
equation, i.e. with the equation 

Sz[X-, z,z)Sz{X-, z,z) — aSz = u{z,z), VA G C (8.31) 

where a = (9^5(0,0) and u = 1 -\- dzS'f^'^. The eikonal equation ()3.26aj) is 
obtained as a reduction of equation (|8.31() taking 5(0, 0) independent on z. 
In fact, the B condition ()8.3)) . producing 5(0,0) = realizes this reduction. 
In what follows we will discuss how one has to choose the 9-data in a way 
to construct the dVN hierarchy directly. We will find constraints which are 
dispersionless analog of the constraints found in |64j . which specify two- 
dimensional Schrodinger equation with real-valued potential. In particular, 
we will see that it is possible to weaken slightly the S-condition ()8.3|) . since 
the value 5(0, 0) is fixed up to a constant by dVN reduction. 

In the following we will focus on solutions of d problem of the form 
5 = 5o + 5, where Sq has polynomial singularities at A = and X = co and 
5 is holomorphic at these points and such that 

lim5(A,A)=0, (8.32) 
lim 5 (A, A) = 5(0, 0, z, z) = -iv{z, z), (8.33) 
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Lemma 1 Let the kernel W in equation (|8.1|) satisfies the assumptions of 
the Vekua's theorem, and let S be its solution. The condition 

S [X,X) = -S {-X,-X) + const (8.34) 

is verified if and only if 

W (A, A, 5a (a, a)) = T^' (-A, -A, -5a (-A, -A)) (8.35) 



Proof. The condition (|8.35|) is a necessary one. Indeed, starting from (|8.1() . 
one obtains 



^ (5, (A, A)) = W (A, A, 5a (a, A)) — (5, (A, A)) (8.36) 

_ (5, (-A, -A)) = W (-A, -A, -5a (-A, -A)) — (5, (-A, -A)) (8.37) 

Exploiting the condition (|8.34|) in (|8.37j) , one gets the equahty (|8.35j) . 
The condition ()8.35|) is sufficient. Indeed, let us introduce the function 

$ (A, A, z, z) = S, (A, A) + 5, (-A, -A) . (8.38) 

Using the equations (|8.32|) and 1)8. 33() . one has 

lim $ (X,X,z,z) = 0, (8.39) 

lim<^ {X,X,z,z) = -2iv^{z,z) . (8.40) 

A — ^0 

Both terms in the right hand side of 1)8.38(1 satisfy the Beltrami equa- 
tion ()8.36|) and ()8.37j) respectively, from which, exploiting the constraint ()8.35|) . 
one concludes that 

^ (A, A, z, z) = W (A, A, 5a (A, A)) ^ (A, A, z, z) . (8.41) 

The function <I> is a solution of the Beltrami equation vanishing at A — > oo. 
So, $ vanishes identically on whole A-plane, so 

f (A,A)=-f (-A.-A). (8.42) 

In particular <^ (0, 0, z, z) = —2ivz {z, z) = ,tliat is Vz {z, z) = 0. Analo- 
gously it is possible to demonstrate that 

§(A,A) = -§(-A,-A). (8.43) 

Hence Vz (z, z) = 0. Equations (jO^ and lead to the relation 

□ 
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Lemma 2 Let S (A, A, z, be a solution of d problem (|8.H) such that (|8.32j) 
and H8.33() are verified. Then the condition 

S[\\)=s{l,^+iv{z,z) (8.44) 

is verified if and only if 

(A, A, 5a (a, ~\)) = -W Q, p -\^S-^ ^ ) . (8.45) 

Proof. The condition (|8.45|) is a necessary one. Let us consider the 
complex conjugation of equation (|8.1() 

dx {S (A, A)) = W (A, A, Sx (A, A)) . (8.46) 

Since 

S (A, A) = 5 + iv{z, z) = S (e, C") + z) (8.47) 

where ^ = A^^ and ^ = A^^, the left hand side of (|8.4(i|l can be written as 
follows 

S>.(S(i.m = §f (f,f) = -^M'(?.?.8£S(tD) = 

that provides us with equation (|8.44|) . 

The condition 1)8. 45(1 is a sufficient one. The 9-equation 1)8. 1|) written in 
terms of the variables ^ = A~^ and ^ = A^^ 

= VF(e,e"5g(e,e)) (8.49) 

is equivalent to 

Using equation (|8.46j) . multiplied by A^, and equation (|8.45|) . one concludes 
that 

dxS {\,\) = dxS {j,'^ . (8.51) 



Integrating (|8.51|) . one gets 



5(A,A) =5(i,i) +^(z,z). (8.52) 
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Let us note that v cannot depend on A since the function S is mero- 
morphic outside the ring. Now, evaluating the equahty H8.52|) at A = and 
A — > c«, one obtains 

V (z, z) = ^(0, 0, z, z) = -5(0, 0, z, z). (8.53) 
So, ?) is a purely imaginary function 

v{z,z) = iv{z,z) (8.54) 

where f is a real valued function. This complete the proof □. 

Note that Vz = ivz{z, z) = —a, in other words, it is the coefficient in front 
of the "magnetic" term in equation (|8.3H) . When it vanishes (i.e. o = 0), 
one has the pure potential equation (|8.31() . that is the eikonal equation. 

Combining together the lemmas and , one gets the following theo- 
rem: Theorem If the d-data W of the d-equation 1)8. 1(1 obey the constraints 

A2T^(A,A,5a(A,A)) =-W^Q,p-A25x(i,i)) , (8.55) 

W (A, A, Sx (A, ~X))=W' (-A, -A, -Sx (-A, -A)) , (8.56) 

then this d-problem provides us with the eikonal equation with real-valued 
refractive index. 

9 Concluding remarks. 

Most results in nonlocal nonlinear optics, except very few interesting cases jSl 
lllj . have been obtained previously by numerical analysis of the nonlocal 
NLS equation (see e.g. ^l^JE]). We have demonstrated the existence of 
integrable regimes for the Cole-Cole media in the geometric optics limit. In 
this case many calculations can be performed analytically. Of course, this 
type of model provides us with certain approximation for the phase of the 
electric field and of intensity. As it happens for the helicoidal wavefronts, 
the model breaks along the lines where the phase is singular. Then, even 
if our approach allows us to analyze exactly the properties of wavefronts 
and the structure of their dislocations, the approximation for intensity is no 
more accurate. The helicoid is a particular meaningful example of this fact. 
Indeed, our geometric optics model gives a blows up instead of a vanishing 
intensity As observed above, the intensity of light beam vanishes due 
to the no more negligible wave corrections along the axis of helicoid, where 
the phase is singular. These vortex type structures are of interest from the 
phenomenological point of view because the are often associated with dark 
solitons in nonlocal medi [HSja. So, one may expect that the singularities 
of the phase may give a clue of possible existence of particular interesting 
structures just like dark solitons. 
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We also have demonstrated that in the Beltrami sector there exist self- 
guided light beams exhibiting nontrivial singular phase structures which 
were not considered before. An accurate description of these wavcfront dis- 
locations could gives useful indications about the physical properties of light 
beams beyond our model. 
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